Ferroelectricity from spin supercurrents in LiCuV04 
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We have studied the magnetic structure of the ferroelectric frustrated spin- 1/2 chain material 
LiCuV04 in applied electric and magnetic fields using polarized neutrons. A symmetry and mean- 
field analysis of the data rules out the presence of static Dzyaloshinskii-Moriya interaction, while 
exchange striction is shown to be negligible by our specific-heat measurements. The experimen- 
tally observed magnetoelectric coupling is in excellent agreement with the predictions of a purely 
electronic mechanism based on spin supercurrents. 

PACS numbers: 77.80.-e, 75.10.Jm, 75.10.Pq, 71.70.Ej 
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Materials where ferroelectricity (FE) appears simulta- 
neously with magnetic order are referred to as type-II or 
magnetic multiferroics. Considerable attention has been 
devoted to their unconventional magnetoelectric (ME) ef- 
fects [l| , namely the control of their electric polarization 
P with an external magnetic field H 0], and likewise at 
H=0 the favoring of a particular magnetic domain by an 
applied electric field E 043- However, the fundamental 
mechanisms responsible for FE in this class of materials 
are under debate, in particular whether FE always results 
from atomic displacements or whether it can be induced 
by a purely electronic mechanism based on spin super- 
currents [8| via the Aharonov-Casher effect, a relativistic 
topological quantum effect for neutral magnetic particles. 
Here, like in the spin Hall effect Q, the electric field di- 
rectly controls the spin current without intermediating 
atomic displacements. 

Phenomenologically, a spin spiral may induce an elec- 
tric polarization P = Aei2 x (Si x S2), where 612 links 
neighboring spin sites Si and S2 and A depends on the 
fundamental mechanism of the ME coupling . Mech- 
anisms that are experimentally established involve sym- 
metry breaking lattice distortions, which are either in- 
duced by symmetric exchange [ll| or b y a syrnmetric 
Dzyaloshinskii-Moriya interaction (DMI) IH . The 
spin supercurrent mechanism Q, on the other hand, op- 
erates in the absence of magnetostriction and static DMI. 
Here frustration of the symmetric exchange interactions 
leads to non-coUinear magnetic order and a spin super- 
current that couples to the electric field via the relativis- 
tic spin-orbit coupling (SOC). The sign and size of A 
is therefore materials specific [1, However, the spin 
supercurrent vanishes if the spin spiral is stabilized by 
DMI Q, and the latter is present in many materials, e.g. 
those studied in Refs. (i-Q. 

In this Rapid Communication, we provide strong ex- 
perimental evidence that the ferroelectricity in the frus- 
trated spin-1/2 chain material LiCuV04 is purely driven 
by the spin supercurrent mechanism and that 
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FIG. 1: (Color online) (a) Unit cell of LiCuVOd where only 
edge-sharing CUO4 chains are shown. Solid (dashed) lines 
represent the ferromagnetic (antiferromagnetic) exchange in- 
teractions, (b) Isolated edge-sharing CUO4 chain with single- 
hole d^2_y2 orbitals and a —im'^ cycloid, (c) Heat capacity at 
the magnetic transition shows absence of hysteresis between 
heating and cooling. 



magnetostriction is negligible. By combining mean-field 
theory and symmetry analysis with polarized neutron 
diffraction measurements of the magnetic structure in 
crossed electric and magnetic fields, we show that the 
static DMI plays no role in driving LiCuV04 to a spin- 
spiral ferroelectric state. 

LiCuV04 crystallizes in the orthorhombic Imma space 
group (#74), where the two Cu^+ ions in the primitive 
unit cell are located at the positions ri = (0,0,0) and 
r2 = (0, 1/2,0) [Fig. [IKa)]. They form spin-1/2 chains 
along the b axis built up from edge-sharing CUO4 units. 
Below T/v = 2.4 K, three-dimensional antiferromagnetic 
(AFM) order sets in with an incommensurate propaga- 
tion vector K=(0, Ky,0) with Kj, = 0.532 [15|, stabilized 
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by frustration of nearest neighbor (Ji) ferromagnetic 
(FM) and next- nearest neighbor (J2) AFM exchange in- 
teractions along the chain [3 12 1- Simultaneously, an 
electric polarization appears along the a axis [isj . In a 
magnetic field along the a axis, a spin-flop transition oc- 
curs at Hf=2.^ T [H, [i^. A concomitant flop of P 
from the a to the c axis is observed for H || a, while no 
electric polarization is reported alon g b i rrespective of the 



magnetic field direction above Hf^ 21 

Heat-capacity measurements show that the magnetic 
order in LiCuV04 develops directly from the param- 
agnetic Imma phase in a single second-order phase 
transition without any hysteresis within the accuracy 
(<5 mK) of the relaxation method of our Quantum De- 
sign calorimeter [Fig. [Ijc)]. This indicates that symme- 
try breaking atomic displacements are negligible, which 
is also confirmed by our single-crystal neutron diffraction 
measurements below Tjv. 

We have determined the zero-field magnetic struc- 
ture of LiCuV04 by reanalyzing the unpolarized neutron 
diffraction data of Ref. [1^. We find that the ordered 
Cu moments lie in the ab plane forming a planar spin 
spiral, which propagates along the b axis. This cycloid 
structure can be described in terms of Fourier compo- 
nents of the magnetic moments at the two Cu positions as 
= (ma, ibimj,', 0) and nij = — mf exp(— ittKj,) with 
nia and m'^ being real and positive. The ± sign corre- 
sponds to two possible chirality domains. 

The population of these two chirality domains can be 
controlled by an electric field, as shown by our spherical 
neutron polarimetry measurements. In this technique, 
the directions of the incoming (P,) and outgoing (P/) 
neutron polarization can be chosen independently, us- 
ing a so-called "cryopad" device This device was 
mounted on the thermal triple-axis spectrometer IN20 
at the Institut Laue-Langevin. Polarized neutrons of 
wave vector 2.662 were produced and analyzed by 
the (111) reflection of Heusler crystals as monochromator 
and analyzer, while second-order (unpolarized) neutrons 
were removed by a PG(002) filter in the scattered beam. 
A single crystal of LiCuV04 of size 1.9 x 0.7 x 5 mm^ 
was aligned in the be plane and cooled down to T~1.5 
K in a vertical electric field Ea — 45 kV/m pointing in 
the —a direction. Scans along (OfcO) were performed for 
magnetic Bragg reflections (0, ±k',±£) with £ = 1,3,5, 
and k' = 1 — Ky = 0.47. Introducing a cartesian coordi- 
nate system with x along Q, y perpendicular to Q in the 
scattering plane, and z perpendicular to the scattering 
plane, the neutron cross-section for P^ along y and P/ 
along ±x is (7.y:r,y2 0ci(|M^Q|2^z[M;5_QxM_LQ]-x). Here, 
M±Q = Q X [M(Q) X Q]/|QP is the component of the 
magnetic structure factor, M(Q) = m^exp(iQ ■ r^), 
that is perpendicular to the wave-vector transfer Q = 
(h, k,£). We note that the chiral term i[M^q xMj^q] is 
either parallel or anti-parallel to Q. 

Figurc[5Ja) shows the resulting (Jyx,yx cross-sections for 
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FIG. 2: (Color online) (a-c) Neutron polarization dependence 
of the magnetic reflections (0, ±0.47, ±1) at T=1.5 K with an 
electric field Ea directed along the —a direction. The lower- 
left panel is a sketch of the Brillouin zone in the be plane. 



the magnetic Bragg peak (0, —k', 1). For the cycloid with 
—im'l sketched in Fig. [TJb), the b component of M(Q) 
is negative (the dashed red arrow shows its projection _L 
Q), which implies that the chiral term i[M^QxMxQ] in 
the cross-section is parallel to Q. Intensity should then 
be observed only in the Uy^ channel, which is the case. At 
Q = (0, k', 1) [Fig. mb)], the direction of the propagation 
vector K is reversed, which in turn reverses the b compo- 
nent of M(Q), since m"'^ = (m^)*. The sign of the chiral 
term is therefore reversed, i.e. it is anti-parallel to Q, and 
intensity is observed in Oyx, still corresponding to the 
same chirality of the cycloid {—im'^). At Q = (0, k', — 1) 
[Fig. [2{c)], the directions of both Q and k. are reversed, 
and the situation of Fig. [2](a) is recovered. Cooling the 
sample in zero electric field leads to equal intensity in 
the two polarization channels, which implies that both 
domains are equally populated. From these observations 
we conclude that an electric field E || a controls the chiral 
domain population and that the ME coupling is nega- 
tive, A<0, since the populated domain has ei2x(SixS2) 
anti-parallel to P (and E). 

To obtain more quantitative information on the mag- 
netic structure for E || a, we modeled all 36 polar- 
ized cross sections (a, /? = ±x, ±y, ±z) of the five 
magnetic Bragg peaks measured using a more general 
= {ma,ni'i^ -f im|,',0), i.e. allowing nih to have a real 
component in addition to the imaginary one considered 
until now. A least-squares fit was performed with m^, 
m'fj, m'^, and d- as free parameters, where d- is the 
domain population of the —im'^ domain. The result- 
ing parameters toJ^ = 0.01(5) m^, m'^ = 0.95(3) ma, and 
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d- =0.99(2) establish that the zero-field magnetic struc- 
ture is a circular cycloid in the ab plane with a single 
chirality domain fully populated. A similar fit of two 
refiections measured with £'=0 gave an equal domain 
population, i.e. d_=d+=0.50(0), which strongly supports 
the conclusion that full control of the cycloid chirality is 
achieved by applying an electric field along a. 

These results are confirmed when a magnetic field H 
is applied close to the c axis, in which case no spin-flop 
transition is expected nor observed. In order to detect 
the chiral term we used the same &c-oriented sample as 
before with a magnetic field in the horizontal scattering 
plane and a vertical electric field E || a. Compared to 
the previously described set-up, we used an unpolarized 
incoming beam from a Si(lll) monochromator but kept 
polarization analysis and a spin fiipper in the scattered 
beam. In this configuration, P/ is parallel or anti-parallel 
to H. The neutron polarization for a sample cooled in a 
positive electric field is shown in Fig. [3ljb), and, follow- 
ing the same reasoning as for Fig. [5J corresponds to a 
cycloid with +im'^. The sample was then warmed up to 
10 K, well above T/v, and cooled down in a reversed elec- 
tric field. This leads to a complete inversion of the chiral 
domain population [Fig. [3l^c)]. Complementary measure- 
ments with H=S.5 T exactly along the c axis corroborate 
these findings [see Fig.[31[d)]. From this, we conclude that 
the ab cycloid is stable upon application of H || c up to 
at least 3.5 T and that an electric field along ±a fully 
controls the chirality of the cycloid. Furthermore, we 
find that the ME coupling is negative, A<0, as in zero 
magnetic field. 

A crucial test of the validity of the spin supercurrent 
model for LiCuV04 is to raise the magnetic field above 
the spin-flop field H^^ , where theory predicts a cycloid 



in the be plane, which according to Ref. jl4| results in a 
positive ME coupling, A> 0. Measurements were per- 
formed on a 4.5 x 3.5 x 1 mm'^ single crystal aligned 
in the ab plane and cooled down with H || a=3.8 T 
and E || +c as sketched in Fig. ^e). A magnetic peak 
was found at Q=(1,k',0), which confirms the stability 
of the propagation vector k for Ha > H^^ . As shown 
in Fig. [Sjf-h), Pf changes sign when E is reversed (un- 
der field-cooled conditions) and essentially vanishes for 
zero electric field. For an electric field E || -|-c (— c) 
of 150 (130) kV/m, we observe a circular cycloid given 
by m^' = (0,m, ±im) with +im {—im) and a domain 
population of rf+=0.87 (d_=0.78) [Fig. [^f-g)], while for 
zero electric field an equal domain population is observed, 
d+=0.55(2) [Fig.EJh)]. We conclude that an electric field 
along the c axis controls the magnetic domain population 
of the spin cycloid observed above the spin-flop field, at 
Ha= 3.8 T. The chirality of the cycloid is such that the 
ME coupling constant A is positive for the 6c-plane struc- 
ture, A> 0. 

Symmetry analysis, which is directly applicable since 
there is only one second order phase transition, leads 
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FIG. 3: (Color online) (a-d) Neutron polarization dependence 
of the magnetic reflection (0,0.47, 1) at r=2.1 K with Ea^^ 
±45 kV/m and ^=1 T along Q as sketched in (a) or H^3.5 T 
along c. (e-g) Similar scans for (1, —0.47, 0) at T=1.9 K with 
Ec ^0 and Ha =3.8 T along a as sketched in (e). The Ec<0 
value had to be reduced with respect to Ec>0 due to electric 
discharges, (h) Same scans for E=0. Identical results are 
obtained with the direction of H reversed. 



to four one-dimensional irreducible representations (IR), 
none of which can describe the observed ab cycloid. 
To proceed, we calculate the mean-field ground-state 
energy, egs{K) = Eij[S*(-'«)J*j(K)Sj (/«)], where 
Jij (k) is a (complex) 6x6 matrix representing the 
Fourier-transformed exchange interactions. Imposing the 
symmetry constraints of the Imma space group on J^j (k) 
leads to cancelation of all off-diagonal terms except terms 
of type yz. This implies that the DMI term of type 
xy (Dz) is zero for every bilinear exchange path, which 
precludes the stabilization of an ab cycloid by DMI. 
A 6c-cycloid cannot be stabilized either by DMI, since 
the yz {Dx) term alternates in sign along the chain. 
If we impose as constraints our experimental findings 
that h+l is odd, TOc=0, and that the magnetic structure 
contains both nia and mf, components, it follows from 
minimizing egs(K) that the exchange matrix is pseudo- 
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tetra gon al for all dominant exchange paths defined in 
Ref. 116|, unless there are pathological coincidences of 
the exchanges. The cycloid is then stabilized in the ab 
plane by a diagonal exchange anisotropy of easy-plane 
type, \Jz\ < \Jx\ = \Jy\- Imposing as maximum spin 
value S=i on both Cu sites, we find that the general 



that the spin supercurrent creates FE polarization in the 
absence of atomic displacements, only by asymmetric or- 



bital occupancy These calculations predict the sign 



Fourier component of site 1, mf = (ttIq, ±[to[, 



,0), 



is reduced to = (m, ±im, 0) with m real. This cor- 
responds to the combination of two irreducible repre- 
sentations with equal amplitude in quadrature and rep- 
resents two degenerate perfectly circular cycloids of ei- 
ther rotation sense, in agreement with our experimental 
findings. From these mean-field considerations, we con- 
clude that the zero-field spin-cycloid of LiCuV04 occurs 
without stabilizing DMI, driven exclusively by frustrated 
anisotropic diagonal (pseudo-tetragonal) exchange. 

A magnetic field applied in the ab plane leads to a 
spin-flopped cycloid or helix in a plane perpendicular to 
the field as soon as the Zeeman energy overcomes the 
anisotropic diagonal exchange. Our finding of a (near) 
single-domain spin-flopped 5c-plane cycloid for Ha > H^^ 
and E II c with the same propagation vector as before and 
a chirality that depends on the sign of the electric flcld 
confirms the purely diagonal exchange scheme, and is 
quite distinct from the complex behavior of other mate- 
rials [l3|. The anisotropy of the diagonal exchange en- 
ergy calculated within mean-field theory using the spin- 
flop field of 2.5 T and the main exchange integrals from 
Ref. [3 is less than ~0.5%. Similar energy consider- 
ations for the spin-flop field put an upper limit on the 
ratio of off-diagonal (yz) and diagonal exchange energy 
to less than 4-10~^. The exchange matrix of LiCuV04 
is hence remarkably isotropic, with diagonal anisotropic 
exchange being the leading (small) perturbation. This 
pseudo-tetragonal spin symmetry of the exchange ma- 
trix has its origin in the almost tetragonal symmetry of 
the CuOg octahedron and the perfectly flat edge-shared 
Cu04 chains in the ab plane of LiCuV04. The latter 
lead to a pseudo-tetragonal exchange matrix even in the 
presence of a small (orthorhombic) distortion of the bond 
angles away from 90° [2^ . 

By combining our experimental results with mean-field 
calculations, we have shown that the formation of a fer- 
roelectric spin-cycloid in LiCuV04 below Tjv is driven by 
purely diagonal frustrated anisotropic exchange without 
stabilizing static DMI or magnetostriction. This rules out 
most of the currently discussed mechanisms for FE except 
the spin supercurrent scenario Q . This purely electronic 
scenario has been studied for both linear Cu-O-Cu clus- 
ters of C n electr ons Q and for edge-sharing CUO4 geom- 



etry 14, 23, 2^1. In the latter case, SOC has been shown 
to generate a small symmetric exchange anisotr opy 12311 
as well as Coulombic ring-exchange processes [231, l25l |. 
which lead to chiral correlations and finally ME coupling 
psj . Tetramer diagonalization and first-principles calcu- 
lations of the SOC on O- and Cu-sites of LiCuV04 show 



change of the ME constant A at the spin-flop transition 
for H II a that we observe experimentally. 

We conclude that the purely electronic scenario of spin 
supercurrents is realized in LiCuV04: a small symmet- 
ric exchange anisotropy is present while magnetostriction 
and static DMI are negligible, and the theoretical predic- 
tions of both the ratio of the electric polarizations Pa/Pc 
and the sign of the ME coupling in all ferroelectric phases 
are found in excellent agreement with experiments. 

We acknowledge fruitful discussions with J. Schweizer, 
G. Nenert, T. Ziman, and M. E. Zhitomirsky and ex- 
perimental assistance from P. Steffcns, P. Chevalier, 
E. Bourgcat-Lami (neutrons), and G. Siegle (specific 
heat). 
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